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Abstract. The enumeration of independent sets of regular graphs is of interest 
in statistical mechanics, as it corresponds to the solution of hard-particle models. 
In 2004, it was conjectured by Fendley et al. that for some rectangular grids, 
with toric boundary conditions, the alternating number of independent sets is 
extremely simple. More precisely, under a coprimality condition on the sides of 
the rectangle, the number of independent sets of even and odd cardinality always 
differ by 1. In physics terms, this means looking at the hard-particle model on 
these grids at activity —1. This conjecture was recently proved by Jonsson. 

Here we produce other families of grid graphs, with open or cylindric boundary 
conditions, for which similar properties hold without any size restriction: the 
number of independent sets of even and odd cardinality always differ by 0, ±1, or, 
in the cylindric case, by some power of 2. 

We show that these results reflect a stronger property of the independence 
complexes of our graphs. We determine the homotopy type of these complexes 
using Forman's discrete Morse theory. We find that these complexes are either 
contractible, or homotopic to a sphere, or, in the cylindric case, to a wedge of 
spheres. 

Finally, we use our enumerative results to determine the spectra of certain 
transfer matrices describing the hard-particle model on our graphs at activity — 1, 
These results parallel certain conjectures of Fendley et al., proved by Jonsson in 
the toric case. 



1. Introduction 

The hard-square model is a famous open problem in statistical mechanics. In this 
model, some of the vertices of an N by N square grid are occupied by a particle, with 
the restriction that two adjacent vertices are never both occupied (Figure[TJ. In graph 
theoretic terms, an admissible configuration of particles is just an independent set of 
the square grid, that is, a set of pairwise non-adjacent vertices. The key question 
is to enumerate these sets by their size, that is, to determine the following partition 
function at activity u: 

zn{u) = y^u^\ 

I 

where the sum runs over all independent sets of the grid. This problem is highly 
unsolved: one does not know how to express Zjv(u), nor even the thermodynamic 
limit of the sequence Z^{u) (that is, the limit of Zn{u) 1 I n ). The most natural 
specialization of Zn(u), obtained for u = 1, counts independent sets of the N x N- 
grid. It is also extremely mysterious: neither the sequence Zjv(1), nor the limit of 
Zn{1) 1 I n (the so-called hard-square constant) are known. We refer the reader to the 
entry A006506 in the On-line Encyclopedia of Integer Sequences for more details [10] . 
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Note that the thermodynamic limit of Zjq{u) is known if one replaces the square grid 
by a triangular one — a tour de force achieved by Baxter in 1980 [T]. 



Figure 1. A hard particle configuration — or an independent set — 
of the 7 x 7- grid. 

In 2004, Fendley, Schoutens and van Eerten [3] published a series of remarkable 
conjectures on the partition function of the hard-square model specialized at u = —1. 
For instance, they observed that for an M x iV-grid, taken with toric boundary 
conditions, the partition function at u = — 1 seemed to be equal to 1 as soon as M 
and N were coprime. They also related this conjecture to a stronger one, dealing with 
the eigenvalues of the associated transfer matrices. These conjectures have recently 
been proved in a sophisticated way by Jonsson [8]. 

One of the aims of this paper is to prove that similar results hold, in greater 
generality, for other subgraphs of the square lattice, like the (tilted) rectangles of 
Figure [2] (they will be defined precisely in Section [3]). For these graphs, we prove 
that the partition function at u = — 1 is always 0, 1 or —1. 



(0,0) 

ft(6,8) 71(6,9) ft(7,8) 7£(7,9) 



Figure 2. The rectangular graphs 1Z(M,N), defined in Section [3l 

We then show that these results actually reflect a stronger property of the inde- 
pendence complex of these graphs. The independent sets of any graph G, ordered by 
inclusion, form a simplicial complex, denoted by S(G). See Figure[3j where this com- 
plex is shown for a 2 x 2-grid G. The (reduced) Euler characteristic of this complex 
is, by definition: 

xg= £ (-i)i'i- 1 . 

Je£(G) 

The quantity |7| — 1 is the dimension of the cell I. The above sum is exactly the 
opposite of the partition function Zq{u) of the hard-particle model on G, evaluated 
at u = — 1. This number, 

Z G (-1) = £ (-l)l'l = -xa, (1) 

will often be called the alternating number of independent sets. The simplicity of 
the Euler characteristic for certain graphs G suggests that the complex S(G) could 
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have a very simple homotopy type (we refer to Munkres [U] for the topological terms 
involved). We prove that this is indeed the case for various subgraphs of the square 
lattice. For instance, for the rectangles of Figure [2j the independence complex is 
always either contractible, or homotopy equivalent to a sphere. Our results rely on 
the construction of certain Morse matchings of the complex S(G). Roughly speaking, 
these matchings are parity reversing involutions on S(G) having certain additional 
interesting properties. 

1 2 

£(3,3) 

4 3 



£(5,5) 



Figure 3. Two graphs, and the independence complex of the top one. 
This complex has reduced Euler characteristic 1, and is homotopic to 
a O-dimensional sphere (two points). The patient reader can check 
that for the "Swiss cross" £(5,5), the reduced Euler characteristic is 
—1. The corresponding complex is homotopic to a 3-dimensional 
sphere. 

Let us now describe the contents of this paper, and compare it to Jonsson's results. 
In Section [21 we first review the needed background of Forman's discrete Morse 
theory [1]. We then describe a general construction of Morse matchings for the 
independence complex of any graph. The matchings we construct are encoded by 
a matching tree. In Sections [3] to we apply this general machinery to determine 
the homotopy type of the independence complex of several subgraphs of the square 
grid: the tilted rectangles of Figures [2] and [5] (Section [3]), the parallelograms of 
Figure [8] (Section [5]), and variations on them (Section 17. All these graphs have 
open boundary conditions (as opposed to the toric boundary conditions of [31 [8]). 
However, in Section [H we identify two sides of the rectangles of Figure [2] to obtain 
rectangles with cylindric boundary conditions. Again, we determine the homotopy 
type of the associated independence complex. Note that Jonsson recently went one 
step further by studying the same tilted rectangles with toric boundary conditions; 
but he was only able to determine the Euler characteristic ([7], Section 7). Our results 
deal with a finer invariant of the complex (the homotopy type) and the proofs are 
simpler, but they also refer to easier graphs; the toric case is at the moment beyond 
reach of our methods. 

Finally, in Section we give background about transfer matrices, and show how 
the results of the previous sections can be used to derive the spectrum, or at least 
part of the spectrum, of several transfer matrices naturally associated to our graphs. 

We conclude the paper with a discussion on possible extensions of our work, with 
the double objective of discovering new subgraphs with a simple alternating number 
of independence set, and addressing the conjectures of Fendley et al. in the cylindric 
case. 



{M} {2,4} 

{1} {2} {3} {4} 
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2. Morse matchings for independence complexes 

The key tool in the proof of our results is the construction of Morse matchings on 
independence complexes. In this section we give the background and explain how we 
use a matching tree as a systematic construction for Morse matchings. 

2.1. Generalities 

Let us first recall that the poset of independent sets of a graph G, ordered by 
inclusion, is a simplicial complex, denoted S(G). 

We regard any finite poset P as a directed graph, by considering the Hasse diagram 
of P with edges pointing down (that is, from large elements to small elements). 

A set M of pairwise disjoint edges of this graph is called a matching of P. This 
matching is perfect if it covers all elements of P. This matching is Morse (or acyclic) 
if the directed graph obtained from P by reversing the direction of the edges in M is 
acyclic. For instance, the matching of the complex of Figure formed of the edges 
(0,{2}), ({3}, {13}), ({4}, {24}) is acyclic. 

Theorem 1 ([5], Theorem 6.3). Let £ be a finite simplicial complex, seen as a poset, 
and M a Morse matching on S, such that the element of £ is matched. For i > 0, 
let rii be the number of unmatched i-dimensional elements o/£. Then there exists a 
CW-complex having 1 + no 0- dimensional cells and n; L i-dimensional cells for i > 1 
that is homotopy equivalent to £. 

Again, we refer to [9] for the topological terms involved. We will only use the 
following immediate corollary. 

Corollary 2. Under the assumptions of Theorem\T\ if all unmatched elements in £ 
have the same dimension i > and there are j of them, £ is homotopy equivalent to 
a wedge of j spheres of dimension i. In particular, if £ is perfectly matched, then it 
is contractible. 

2.2. Matching trees 

Let us now describe the general principle we use to construct Morse matchings for 
a complex £(G). Let V denote the vertex set of G. The most naive way to define a 
matching of £ = £(G) is probably the following. Take a vertex p £ V, and denote 
by N(p) the set of its neighbours. Define 

A = {/ e £ :/nJV(p) = 0}. 

The set of pairs (I, I U {p}), for / G A and p £ I, forms a perfect matching of A, 
and hence a matching of £. We call p the pivot of this matching. The unmatched 
elements of £ are those containing at least one element of N(p). There may be 
many unmatched elements, but we can now choose another pivot p' to match some 
elements of S \ A, and repeat this operation as long as we can. Of course, the 
resulting matching will depend on the successive choices of pivots. 

This rather naive idea is the leading thread in the construction of our Morse 
matchings of £. In some occasions, we will have to split the set of yet unmatched 
elements into two subsets, and choose a different pivot for each of them. This explains 
why our matching procedure will be encoded by a branching structure, namely a plane 
rooted tree, called a matching tree of £. The nodes of this tree represent sets of yet 
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unmatched elements. Some nodes are reduced to the empty set, and all the others 
are subsets of E of the form 

E(A, B) = {/GS:iC/andBn/ = 0}, 

where 

A n B = and N(A) := U aeA N{a) C B. (2) 

We say that the vertices of A U B are the prescribed vertices of T,(A,B). The root 
of the tree is S(0,0) = S, the set of all independent sets of G. The sets A and B 
will increase along branches, making the sets £(A, B) of unmatched elements smaller 
and smaller. The leaves of the tree will have cardinality or 1, and will contain the 
elements that are left unmatched at the end of the procedure. 

Consider a node of the tree of the form E(A, B). How can we match its elements? 
If the node has cardinality 1, that is, if A U B = V, then we are stuck, as there is no 
non-trivial matching of a graph reduced to one vertex. If A U B C V, we may match 
some of the elements of T,(A,B). Pick a vertex p in V' := V \ (A U B). Because 
of J2J, the neighbours of p are either in B, or in V'. This makes p a good tentative 
pivot. If we actually use p as a pivot to match elements of £(A, B), we will be left 
with the following set of unmatched elements: 

C/ = {/GS:AC/,/nS = 0,/n N(p) ^ 0}. 

If p has no neighbour in V' , the above set is empty, and we have perfectly matched 
T*(A, B). If p has exactly one neighbour in V', say v, then U = T>(A U {v }, B U N(v)). 
However, if p has at least two neighbours in V' , say v and v' ', the set U is not of the 
form T,(A',B'). Indeed, some of the unmatched sets / contain v, some others don't, 
but then they have to contain v' . This puts us into trouble, as we want to handle 
only unmatched sets of the form T,(A', B'). We circumvent this difficulty by splitting 
the original set T,(A, B) into two disjoint subsets of the form T*(A' , B'), that differ 
by the status of, say, the vertex v. More precisely, we write: 

S(A,B) = E(A,Bu{v}) W E(AU{v},BUN(v)), 

and then study separately each subset. 

The above discussion justifies the following construction of the children of a node. 
If this node is the empty set (no unmatched elements), we declare it a leaf. Otherwise, 
the node is of the form B). If A U B = V, the node has cardinality 1, and we 
also declare it a leaf. We are left with nodes of the form S(A, B), with A U B C V . 
Choose a vertex p (the tentative pivot) in V' = V \ (A U B), and proceed as follows: 

• If p has at most one neighbour in V, define A(A, B,p) to be the subset of 
T,(A,B) formed of sets that do not intersect N(p): 

A(A, B,p) = {I G S : A C / and B n / = / n N(p) = 0}. 

Let M(A,B,p) be the perfect matching of A(A,B,p) obtained by using p 
as a pivot. Give to the node T,(A,B) a unique child, namely the set U = 
T,(A, B) \ A(A,B,p) of unmatched elements. This set is empty if p has no 
neighbour in V . In this case, we say that p is a free vertex of X(A, B). lip 
has exactly one neighbour v in V, then U = S(iU{u},BUJV(t))). Index 
the new edge by the pivot p. We say that the 3-tuple (A, B,p) is a matching 
site of the tree. 
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• Otherwise, let us choose one neighbour v of p in V . The node T,(A, B) has 
two children, only differing by the status of v. More precisely, the left child 
is E(A, B U {v}) and the right child is Y,(A U {v}, B U N(v)). The union of 
these two sets is T,(A,B). Label the two new edges by the splitting vertex v. 
We say that (A, B, v) is a splitting site of the tree. 

Observe that the new nodes satisfy Conditions ([2]), unless they are empty 

Given a sequence of choices of tentative pivots and splitting vertices, we obtain a 
matching M of £ by taking the union of all partial matchings M(A, B,p) performed 
at the matching sites of the tree. The unmatched elements are those sitting at the 
leaves of the tree. 

The above construction is rather natural, and we invite the reader to practice 
with the example given in Figure HI In this figure, every (non-empty) node is de- 
scribed by the vertices of A (in black) and B (in white). At the matching site 
(0,0, 1), the elements and {1} (among others) are matched. At the matching site 
({2}, {1, 3, 4, 6}, 5), the elements {2} and {2,5} get matched, among others. At the 
end of the matching procedure, the independent set {2, 7} is the only unmatched 
element of S(G). 

D matching site 



splitting site 
(tentative pivot = 5) 

6 



matching site 




matching site 
(free pivot at 7) 



Figure 4. A subgraph G of the square grid, and one of its matching trees. 

We now aim at showing that the matchings obtained with the above procedure 
are in fact Morse. The following lemma gathers some properties of this construction. 

Lemma 3. Every matching tree satisfies the following properties: 

(1) For every matching site {A, B,p), the matching M(A, B,p) is a Morse match- 
ing of A(A, B,p) (still ordered by inclusion). 

(2) Let (A, B,p) be a matching site with a non-empty childT l (A U {v}, B U N(v)). 
Let I G A(A, B,p) and JG Z(AU {v}, B U N(v)). Then J £ I. 

(3) Let (A,B,v) be a splitting site, I G S(^4,i?U{f}) and J G 
Z(AU{v},BUN(v)). Then J £ I. 

Proof. (1) Consider the Hasse diagram of the poset A(A,B,p) and its directed ver- 
sion, with all edges pointing down. Now, reverse the edges of M(A,B,p). The up 
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edges join two elements of the form / \ {p}, I, so they correspond to adding the 
vertex p. The down edges correspond to deleting a vertex different from p. Clearly 
there cannot be a directed cycle in A(A,B,p). 

(2) The set J contains v, a neighbour of the pivot p, while none of the matched 
sets / of A(A, B,p) contains v. 

(3) Here again, J contains v, but I doesn't. a 

The following easy lemma appears as Lemma 4.3 in Jonsson's thesis [6]. 

Lemma 4. Let V be a finite set and A = Ai ttl A2 a collection of subsets of V , 
ordered by inclusion. Assume that if a E Ai and r E A2 then t £ a. Then the union 
of two acyclic matchings on Ai and A2 respectively is an acyclic matching on A. 

We can now establish the main result of this section. 

Proposition 5. For any graph G and any matching tree ofG, the matching o/S(G) 
obtained by taking the union of all partial matchings M(A,B,p) performed at the 
matching sites is Morse. 

Proof. We will prove by backwards induction, from the leaves to the root, the fol- 
lowing property: 

For every node r of the matching tree, the union of the partial match- 
ings performed at the descendants of r (including r itself) is a Morse 
matching of r. We denote this matching UM(r) (for Union of Match- 
ings) . 

The leaves of the tree are either empty sets or singletons, endowed with the empty 
matching, which is Morse. This supplies the induction base. Consider now a non-leaf 
node of the tree, of the form r = £(A, B). 

Assume (A,B,p) is a matching site. By Lemma Ell, M(A, B,p) is Morse. If the 
(unique) child of r is empty, then UM(r) = M(A, B,p) and we are done. If this 
child is t' = T,(A U {v}, B U N(v)), the induction hypothesis tells us that UM(t') is 
Morse. By Lemma[3l2, we can apply LemmalUwith Ai = A(A,B,p), A2 = t' and 
A = t, where the partial matchings on Ai and A2 are respectively M(A, B,p) and 
UM(t'). This shows that UM(t) is Morse. 

Assume (A,B,v) is a splitting site. By induction hypothesis we already have 
Morse matchings on both children of r, namely Ai = T,(A, B U {v}) and A2 = 
S(j4 U {v}, B U N(v)). Again, Lemma 03 3 allows us to apply Lemma HJ and this 
shows that the union UM(t) of CTM(Ai) and UM(A 2 ) is Morse. 

This completes the induction. The case where r is the root of the tree gives the 
proposition. a 

3. Rectangles with open boundary conditions 

In what follows, we consider I? as an infinite graph, with edges joining vertices 
at distance 1 from each other. For M,N > 1, let 1Z(M,N) be the subgraph of Z 2 
induced by the points (x, y) satisfying 

y<x<y + M — 1 and — y < x < —y + N — 1. 

Examples are shown on Figure El Note that 7L(M, N) contains [^] vertices. Other 
rectangular shapes arise when we look at the subgraph TZ(M, N) of Z 2 induced by 
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the points (x, y) such that 

y < x <y + M - 1 and -y + l<x<-y + N 

(see Figure [5]). More precisely, the graphs TZ(2M + 1, 2N + 1) are not isomorphic to 
any of the TZ(K,L) (the other graphs TZ(M,N) are isomorphic to an 1Z(K,L)). 



(i,o)- 



iz(6, 



6(7,7) 



Figure 5. The graphs K(M,N). 

We study the independence complexes of the rectangles TZ(M,N) and 1Z(M,N). 
Recall the general connection ([T|) between the (reduced) Euler characteristic of these 
complexes and the alternating number of independent sets. We use below the nota- 
tion Zfc(M,N) rather than Znng,N)- 

Theorem 6. Let M,N > 1. Denote m = [M/3] and n = \N/3\. 

• If M = 3 1 or N = 3 I, then E(U(M,N)) is contractible and Z n (M,N) = 0. 

• Otherwise, E(7£(M, N)) is homotopy equivalent to a sphere of dimension 
mn - 1, and Z n {M, N) = {-l) mn . 

The above holds also when replacing TZ(M,N) by 7Z(M,N). 

Remark. We will show in Corollary [TO] that, for N = 3 1 and M > 2 1+ W 2 1, the 
alternating number Zfc(M, N;C, D) of independent sets of 1Z(M,N) having border 
conditions C and D on the two extreme diagonals of slope 1 is actually for all 
configurations C and D. This will indirectly follow from the study of tilted rectangles 
with cylindric boundary conditions performed in Section [H 

Proof. We study the graphs 1Z(M, N) and TZ(M, N) together. We construct Morse 
matchings of the independence complexes of these graphs by following the general 
principles of Section [2l We need to specify our choice of tentative pivots and splitting 
vertices. Our objective is to minimize the combinatorial explosion of cases, that is, 
the number of splitting sites. Consider a node T,(A,B). 

(1) In general, the tentative pivot p = (io, Jo) is chosen in V' = V \ (A U B) so 
as to minimize the pair (i + j, i) for the lexicographic order. That is, p lies 
as high as possible on the leftmost diagonal of slope —1. However, if there 
is at least one free vertex (that is, a vertex of V = V \ (A U B) having no 
neighbour in V) on the next diagonal i + j = io + jo + 1, then we choose one 
of them as the pivot: the only child of £(^4, B) is then the empty set. 

(2) If the tentative pivot has several neighbours in V' = V \ (A U B), then it has 
exactly two neighbours in V', namely its North and East neighbours. Indeed, 
the other two neighbours are smaller than p for the lexicographic order, and 
thus belong to A U B. Take v, the splitting vertex, to be the East neighbour 
of p. 



ON THE INDEPENDENCE COMPLEX OF SQUARE GRIDS 9 

These conventions are illustrated by an example in Figure El As before, the elements 
of A and B are indicated by • and o respectively. Rather than labeling the vertices 
of the graph and the edges of the tree, we have indicated the tentative pivots by *, 
and the splitting vertices by A. 



* 




K(M, N - 3) 



Figure 6. The top of the tree describing the matching of £(7£(8, 6)). 

We are going to prove by induction on ./V and M the following properties, valid 
both for the graphs K(M,N) and K(M,N): 

(A) if M =3 1 or N =3 1, then there is no unmatched element in S, 

(B) otherwise, there is a unique unmatched element, of cardinality mn (and thus 
dimension mn — 1). 

By Proposition El the matchings we obtain are Morse. Hence Theorem [6] follows from 
Properties (A) and (B) using Corollary [2l Let us now prove these properties. We 
first study small values of N. 

(1) If N = 1, the graph is formed of isolated points. Hence the matching 
M(0, 0,p) performed at the root of the tree is a perfect matching. Prop- 
erty (A) follows. 

(2) For iV = 2, we leave it to the reader to check (A) and (B) for M = 1,2,3, 
both for the graphs 7Z(M,N) and TZ(M,N). We then proceed by induction 
on M, for M > 4. 
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(a) For the graph 7Z(M, 2), the root is a matching site, and its unique child 
is E(A,B) with A = {(1,0)} and B = {(0, 0), (1, -1)}. The graph 
obtained by deleting the vertices of AuB is (a translate of) 1Z{M — 3, 2). 
This shows that every unmatched element of £(7£(M, 2)) is obtained by 
adding the vertex (1,0) to (a translate of) an unmatched element of 
E(K(M - 3, 2)). The result then follows by induction on M. 

(b) For the graph 7£(M, 2), the root is a splitting site. Its right child is 
perfectly matched using the free pivot (1, 1), and has the empty set as 
its unique child. The left child of the root is partially matched using 
the pivot (1,0). The unmatched elements are those of £(^4,.B), with 
A = {(1,1)} and B = {(1, 0), (2, 0)}. The graph obtained by deleting 
the vertices of A U B is (a translate of) 7Z(M — 3,2). This shows that 
every unmatched element of T,(7Z(M, 2)) is obtained by adding the vertex 
(1, 1) to (a translate of) an unmatched element of Y,{1Z{M — 3,2)). The 
result then follows by induction on M . 

The study of the case N = 2 leads us to introduce a notation that will be useful 
in our forthcoming inductions. 

Notation 1. Let V be a subset of vertices of the square grid, and assume that there 
exist i,j £ Z such that V is the disjoint union V = V\ l±l (V2 + (where V2 + 

= {v + (i,j) : v G V2}). Let Xi and X2 be two collections of sets on the ground 
sets V\ and V2, respectively. Then Y := {Ii ttl (I2 + (i, j)) : I\ G X\,Ii £ X2} is a 
collection of sets on the ground set V . We use the notation Y = X\ * X2 to denote 
that the elements ofY are formed by the concatenation of an element of X\ with (the 
translate of) an element of X2. In particular, \Y\ = | |-X"2 1 - 

For instance, if U(M,N) (resp. U(M,N)) denotes the set of unmatched elements 
in T,(TZ(M,N)) (resp. Y,(1Z(M, N))), the above observations can be summarized by 

U(M, 2) ^ 17(3, 2) * U(M - 3, 2) and U(M, 2) ^ U(3, 2) * U(M - 3, 2). 

We now return to our induction. 

(3) The case N = 3 is very similar to the case N = 2. One first checks that the 
result holds for M = 1,2, 3, both for the graphs K(M, N) and H(M, N). For 
M > 4, the result is proved by induction on M, after observing that 

U(M, 3) ^ U(3, 3) * U(M - 3, 3) and U(M, 3) ^ U(3, 3) * U(M - 3, 3). 

The following three observations will be useful in the rest of the proof. Firstly, 
the tentative pivot is never taken in the third diagonal of slope —1. Secondly, 
when it is taken in the second diagonal, it is a free pivot. Finally, when the 
matching tree has a non-empty leaf T*(A,B) (that is, when M ^3 1), all 
vertices on the third diagonal belong to B. 

(4) Now, let iV > 4. The key observation is that the top of the matching tree 
coincides, as far as the prescribed vertices, pivots and splitting sites are con- 
cerned, with the matching tree obtained for N = 3. This is illustrated in 
Figure and holds as long as the pivots are taken in the first two diago- 
nals. Once these pivots have been exhausted, we are left with (at most) one 
non-empty unmatched set T,(A,B), whose prescribed vertices are those of 
U(M, 3) (if we work with K(M, N)), or U(M, 3) (if we work with TZ(M, N)). 
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Moreover, the tree rooted at the vertex T,(A, B) is isomorphic to the matching 
tree of TZ(M, N - 3) (resp. 1Z(M, N - 3)). This leads to 

U(M, N) ^ U(M, 3) * U(M, N - 3), 

U(M, N) ^ U(M, 3) * U{M, N - 3) 

when for M =3 1 both U(M, 3) and U(M, 3) are empty. Properties (^4) and 
(.B) easily follow. 

■ 

4. Rectangles with cylindric boundary conditions 

We now study a "cylindric" version of the graphs 1Z{M, N), obtained by wrapping 
these graphs on a cylinder. For M, N > and M even, we consider the graph 
TZ C (M, N) obtained from 7L(M + 1,N) by identifying the vertices and (M/2 + 
i, -M/2 + i), for < i < [^f 1 ] . Observe that the rectangles K(M, N) of Figure (5J 
when wrapped in a natural way around a cylinder, yield the same family of graphs. 
We denote by Z^(M, N) the alternating number of independent sets on the graph 
K C {M,N). 

Theorem 7. Let M,N > 1, with M even. Denote m = L^p-J and n = \N/3\. 

• IfN =3 I, then Y,(K C (M , N)) is contractible and Zfc(M,N) = 0. 

• Otherwise, 

— If M =3 0, then £(7£ C (M, N)) is homotopy equivalent to a wedge of 2 n 
spheres of dimension mn — 1, and Z%(M, N) = 2". 

— If M =3 1 or 2 then Y,(1Z C {M,N)) is homotopy equivalent to a single 
sphere of dimension mn — 1, and Z^(M, N) = (— l) n . 

Proof. We define a matching of E = £(7£ C (M, N)) by adopting the same choice of 
tentative pivots and splitting vertices as in the proof of Theorem [6) 
We are going to prove by induction on N the following properties: 

(A) if N =3 1, then there is no unmatched element in S, 

(B) otherwise, 

(Bi) If M =3 0, there are 2 n unmatched elements, each of cardinality mn. 

(B2) If M =3 1 or 2 there is a unique unmatched element, of cardinality mn. 
The theorem then follows from Proposition [5] and Corollary [2l Properties (A) and 
(B) are proved by induction on N. 

(1) When N = 1, the graph is formed of isolated vertices, and the first pivot 
already matches T,(TZ C (M, N)) perfectly. 

(2) When N = 2, the graph is a ring of 2M vertices. Depending on the value 
of M modulo 3, one finds two (if M =3 0) or one (if M =3 1,2) unmatched 
cell(s) of cardinality m. These results are easily obtained by constructing the 
whole matching tree of the graph. We omit the details which are very similar 
to the proof of Theorem [6j 

(3) The case N = 3 is almost identical to N = 2, since the pivot never appears 
in the third diagonal (Figure [7}. As in the case of rectangles with open 
boundary conditions, for each non-empty leaf T,(A,B), all vertices on the 
rightmost diagonal belong to B. There are at most two such leaves (exactly 
two when M =3 0). 
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1l c {M,N-3) 

Figure 7. The top of the matching trees of £(7£ c (6, 5)) and 
£(7£ c (8, 5)). This figure illustrates what happens for M =3 and 
M =3 2. We leave it to the reader to practice with the case M =3 1. 



(4) For N > 4, the top of the tree coincides again with the matching tree of 
1Z C {M, 3). Once the pivots of the first two diagonals have been exhausted, 
the prescribed vertices are exactly those of the first three diagonals. Denoting 
by U C (M, N) the set of unmatched elements of S(^ C (M, N)), it follows that: 

U C (M, N)) ^ U C (M, 3) * U C (M, N - 3). 

Properties (A) and (B) easily follow by induction on N. a 

5. Parallelograms with open boundary conditions 

For K, N > 1, consider now the subgraph V(K,N) of the square grid induced by 
the vertices (x, y) satisfying 

<y < K - 1 and -y<x<-y + N -\. 

An example is shown on Figure El We denote by Zp(K,N) the alternating number 
of independent sets on the graph V(K,N). 

Theorem 8. Let K,N > 1. Denote m = \2K/3]. 
• If K =3 1, then 
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(0,0) I I I I I 

Figure 8. The parallelogram graph "P(4, 8). 

- if N = 3 1 then Y,(V(K,N)) is contractible and Z P (K,N) = 0, 

— otherwise, S("P(J^, N)) is homotopy equivalent to a sphere of dimension 
mn - 1, with n = \N/3], and Z V (K,N) = (-l) n . 

• IfK= 3 2, write N = 2qK + r, with < r < 2K - 1. 

- Ifr = 3 1,2, then E(V(K,N)) is contractible and Z V (K,N) = 0, 

— otherwise ^(^(K, N)) is homotopy equivalent to a sphere of dimension 
mn-l with n = ■ f], and Z V (K,N) = 1. 

• IfK= 3 0, write N = 2q(K + 1) + r with < r < 2K + 1. 

- // r = 3 with r > 1, or r = 3 1 with r < 2K , then E(V(K,N)) is 
contractible and Z-p(K,N) = 0, 

— otherwise, Y,(V(K, N)) is homotopy equivalent to a sphere of dimension 
mn — 1 where n = 1 ^+2 ' lFl> an ^ Z-p{K-,N) = 1. 

Remark. We prove in Corollary [12] that for JV e 3 1 and K large enough, the al- 
ternating number Z-p(K, N; C, D) of independent sets on the parallelogram V(K, N) 
having prescribed conditions C and D on the top and bottom row is 0, for all con- 
figurations C and D. This is not in contradiction with the above theorem: for K 
large enough, if K ^3 1, the quotient q appearing in the theorem is simply 0, so 
that the condition iV =3 1 boils down to r =3 1, with r small compared to K, and 
Z V (K,N) = 0. 

Proof. We construct a Morse matching of the graph V(K, N) by applying the general 
method of Section (2) We then show that for all values of K and N, the matching thus 
obtained has at most one unmatched cell, of cardinality mn. As before, Corollary [2] 
completes the proof. 

The results in the case K =3 1 are reminiscent of what we obtained for rectangles 
(Theorem E]). The rule for choosing tentative pivots and splitting vertices is the same 
as before, and the proof follows the same principles. We do not repeat the argument. 

The other two congruence classes of K are more complicated, and require a differ- 
ent pivot choice. Let us begin with the case K =3 2. First, we partition the set of 
vertices of V(K, 00) into triangular subsets T±, T-j, . . . defined by: 

T 2e+1 = {(x, y): < y < K, 21K <x + y, x-y < 21K + 1 }, 

Tu = {(x, y): 0<y < K, 2(1 - l)K + 1 < x - y, x + y < 21K }. 

These triangles are shown in Figured] for K = 5. 

We now describe the pivot order. The pivots are first taken in 7"i, then in T2 
(once all vertices of T\ are prescribed), and so on. The pivot order for T\ is very 
similar to what we have done previously. That is, we follow the diagonals of slope 
— 1 from upper left to lower right but with the restriction that we must stay within 
the triangle. For T2 the pivot order is what we get when turning T\ upside down. 
That is, we follow the diagonals of slope 1 from bottom left to top right. 
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Figure 9. The configuration of prescribed vertices when the pivots 
of T\ and T 2 have been exhausted, for K = 5. This configuration oc- 
curs in the non-contractible branch of the matching tree. The pivots, 
indicated by * and numbered by pivot order, and the splitting vertices 
A, are empty. Each time a partial matching is performed, some new 
prescribed vertices appear: they are joined by thick line. 



Denote by U P (K,N) the set of unmatched elements of Y>{V{K,N)). We observe 
(Figure EJ) that, once all pivots in T\ and T 2 have been exhausted, only one node 
T*(A, B) of the matching tree is non-empty. The prescribed vertices of this node are 
those of A U B = 71 U T 2 . This gives, for N > 2K: 

U P (K, AO) = U P (K, 2K) * U P (K, N-2K), 

which is the key for our induction on N. By convention, V(K, 0) is the empty 
graph, and its unique independent set is the empty set, so that \U P (K, 0)| = 1. 
We also note that the unique unmatched element of T,(V(K,2K)) has cardinality 
m (m - 1) = (2K + 2)(2K - l)/9. 

Upon iterating the above identity, we obtain, if N = 2qK + r, 

U P (K, AO) U P {K, 2qK) * U P {K, r), 

where the only unmatched cell of Yj{V{K,2qK)) has qm(m — 1) vertices. It thus 
remains to describe what our matching rule produces for the graphs V(K,r), for 
< r < 2K — 1. The following properties are easily observed on the example of 
Figure [9l 

(1) If r =3 1, the rightmost vertex in the top row of V(K, r), which belongs to the 
triangle T\, becomes a free pivot at some stage of the matching procedure, so 
that \U p (K,r)\ = 0. 

(2) If r =3 2, the rightmost vertex in the top row of V(K, r) n 7^ becomes a free 
pivot at some stage, so that again, \U P (K, r)| = 0. 

(3) Finally, if r =3 0, we obtain a unique unmatched cell, of cardinality mr/3. 

Putting together our recursion and the above results for V(K,r), we find that the 
only non-contractible cases are when r =3 0. In this case, there is only one unmatched 
element in T,(V(K, N)), of cardinality qm(m — 1) + mr/3. The result follows for the 
case K =3 2. 
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Figure 10. The configuration of prescribed vertices when the pivots 
of T\ and T2 have been exhausted, for K = 6. 



Let us now adapt this to the final case K =3 0. The triangles that we used to 
define the pivot rule now become trapezoids 71,72, . . ., defined by: 

%i+x = {(x,y): 0<y<K, 2l(K + 1) < x + y, x - y < 2l(K + 1) + 2 }, 

%t = {(x,y): 0<y<K, 2(1 - 1)(K + 1) + 2 < x - y, x + y < 2£(K + 1) }. 

These trapezoids are shown in Figure \W\ for K = 6. The pivots are then chosen 
using the same rule as in the case K =3 2. Again, one observes that when the 
pivots of T\ and T2 have been exhausted, only one vertex T,(A, B) of the matching 
tree is non-empty, and its prescribed vertices are those of T\ U 7^. This gives, for 
N > 2K + 2: 

U P (K, iV)) = U P (K, 2K + 2) * U P (K, N -2K - 2). 

Moreover, S("P(AT, 2K + 2)) has a unique unmatched cell, of cardinality m(m + 1). 
Let us write A = 2q(K + 1) + r, with < r < 2K + 1. Iterating the above identity 
gives 

U P (K, A)) IP(K, 2q{K + 1)) * U p (K,r), 

where the only unmatched cell of T,(V(K, 2q(K + 1))) has cardinality qm(m + 1). 

It remains to describe what our matching rule produces for the graphs V(K,r), 
for < r < 2K + 1. We refer again to Figure fTTJl 

(1) If r =3 1, with r < 2K, the rightmost vertex in the top row of V(K,r) 
becomes a free pivot at some stage of the procedure, so that \U P (K, r)\ = 0. 

(2) If r = 2K + 1, there is a unique unmatched cell, with cardinality m(m + 1) 
(it coincides with the unmatched cell obtained for V(K, 2K + 2)). 

(3) If r =3 with r > 0, the rightmost vertex in the top row of V(K, r) n T2 
becomes a free pivot at some stage of the matching procedure, so that again, 
\U p (K,r)\ = 0. 

(4) If r = 0, we have the empty graph, with the empty set as unique (and 
unmatched) independent set. 

(5) Finally, if r =3 2, we obtain a unique unmatched cell, of cardinality m(r + 
l)/3. 

Putting together our recursion and the above results for V(K, r), we find that the only 
non-contractible cases are when r = 0, r = 2K + 1 and r =3 2. In these cases, there 
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is only one unmatched element in T,(V(K, N)), of cardinality qm(m + 1) + m\r/3\. 
The result follows for the case K =3 0. a 

Remark. The parallelogram V(K, N) gives rise to two distinct families of shapes 
with cylindric boundary conditions: 

• Gluing the two diagonal borders of V(K,N + 1) by identifying the points 
(-i, i) and (-i + N, i) for < i < K - 1 gives the "usual" cylinder Z/iVZ x 
{0, 1,. . . ,K — 1}. It is conjectured in [8| that for odd N the corresponding 
alternating number of independent sets is 1, except for the case N =g 3, K =3 
1 when it is conjectured to be —2. 

• Gluing the two horizontal borders of V(K + 1,N) by identifying the points 
(i, 0) and (i — K, K) for < i < N — 1 gives 1Z C (2K, N), the cylindric version 
of the rectangle which we studied in Section [U 

6. Transfer matrices 

6.1. Generalities 

We develop here a general (and very classical) transfer matrix framework which 
we will instantiate later to the enumeration of independent sets on various subgraphs 
of the square lattice. See [HI Ch. 4] for generalities on transfer matrices. 

Let r > 1, and let S be a collection of subsets in [1, r] := {1,2,..., r}, of cardinality 
d. Let T be a square matrix of size d, with complex coefficients, whose rows and 
columns are indexed by the elements of S. The entry of T lying in row C and 
column D is denoted T(C,D). Call configuration of length n any sequence I = 
(Cq, Ci, . . . , C n ) of subsets of S. We say that Co and C n are the borders of /. The 
weight of I is 

n-l 

w(i) = []T(a,c w ), 

A cyclic configuration of length n is a configuration / = (Co, C%, . . . , C n _i, C n ) such 
that Co = C n . Observe that for all C G S, I = (C) is a cyclic configuration of length 
0, so that there are exactly d such configurations. 

Let t be an indeterminate. It is well-known, and easy to prove, that the length 
generating function of configurations with prescribed borders C and D, weighted as 
above, is 

Gc, D {t)--=Y. E t n w(I) = (l-tT)- 1 (C,D). (3) 

n>0/=(C,Ci,...,C„_i,D) 

From this, one derives that the length generating function of cyclic configurations is 
the trace of (l-fT) -1 : 

G c (t) :=^G c , C (t)=tr(l-tT)- 1 . 

c 

In what follows, we will be interested in deriving eigenvalues of the transfer matrix 
T from the generating functions Gc,D(t) and G c (t). This is motivated by the con- 
jectures of Fendley et al. on the spectra of various transfer matrices related to the 
enumeration of independent sets [3]. One first observation is that finding the whole 
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spectrum (Ai, . . . , A^) of the transfer matrix is equivalent to finding the generating 
function of cyclic configurations: indeed, 

d 

G c (t) = tr(l - tT)' 1 = £Vtr(TT) = £>(A? + ■ ■ ■ + X n d ) = E I^X? ( 4 ) 

n>0 n>0 i=l * 

From PJ, (|3J), and the classical formula giving the inverse of a matrix in terms of its 
determinant and its comatrix, one concludes that 

• For every pair C, D, the reciprocals of the poles of Gc are eigenvalues of 
T, 

• Conversely, the set of non-zero eigenvalues of T coincides with the set of 
reciprocals of poles of the series Gc,c(t), for C running over S. 

In other words: counting configurations with cyclic boundary conditions gives the 
whole spectrum; at least partial information can be derived from the enumeration of 
configurations with open boundary conditions. 

This general framework will be specialized below to the case where the sets C and 
D describe hard-particle configurations (a.k.a. independent sets) on certain layers 
on the square lattice. The weights T(C, D) will be designed in such a way the weight 
of a configuration I is if I is not an independent set, and (— l)' 7 ' otherwise. We 
have schematized in Figure [TT] the various transfer matrix we consider. They will be 
defined precisely in the text. 




Figure 11. How the various transfer matrices act. 



6.2. TWO COMPLETE SPECTRA 

In this section, we combine the above generalities with the results obtained in 
Section 0] for the alternating number of independent sets on the cylinders TZ C (M, N). 
As these cylinders can be generated from two types of transfer matrices (called Mat 
and Lat in Figure [TT1) . we obtain the complete spectra of these two families of transfer 
matrices. All the non-zero eigenvalues are found to be roots of unity. 

We begin with a transfer matrix Pat that describes the independent sets of 
the 2-diagonal graph TZ(2,N). The rows of Pat are indexed by subsets C of 
{1,2,..., [iV/2]}, its columns are indexed by subsets D of {1, 2, ... , [N/2\} : and 

*s{C,D) = lf MDl xcnD = cn(D + i) = 9 t (5) 

otherwise. 
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The notation D + 1 means {i + 1 : i € D}. Observe that Pat is not a square matrix 
if N is o However, if P^r denotes the transpose of Pat, then M.^ := PjvPtv 
is the (square) transfer matrix corresponding to the graph 1Z(3,N): for C,D C 

{i,...,rjv/2i}, 

R N (C,D)= W N (C,E)F N (E,D) = iW + W ^ (-1)'^ 

EC{l,...,[N/2\} E:I=(C,E, D)ind.set 

is, up to the factor il^l+l- !, the alternating number of independent sets I of 7Z(N, 3), 
with top and bottom borders C and D respectively. The coefficients of this matrix 
have actually a simpler expression. Indeed, the configuration (C, E, D) is an inde- 
pendent set if and only if E is in the complement of CUDU (C - 1) U (D - 1). 
Hence the sum of terms (— 1)\ E \ is unless this complement is empty. That is, for 
C,DQ{l,2,...,\N/2}}, 

R N ( C , D) = { i|C|+|D ' if f 1 ' 2 > ■ • ■ > W2J } = C U D U (C - 1) U (D - 1), (g) 
\ otherwise. 

Prom one derives that the entry (C, C) in the product = (PatPaO^ is the al- 
ternating number of independent sets on the cylinder 1Z c (2k, N) studied in SectionlH 
with border condition C on the first diagonal. By Section I6.ll these numbers are 
related to the spectrum (Ai, . . . , \d) of the transfer matrix Mat. More precisely, (J4j) 
gives: 

d 

G c (t) = d + Y, ZWk, N)t k = 77TT! (7) 

fe>l i=l 1 

where the numbers Z^(2k, N) are given in Theorem [7] and d = 2^ N ^ 2 \ 

Theorem 9. Let N > 1. The transfer matrix Mat defined by JB)) has size 2^ N / 2 \ 
If N =3 1, then Mat is nilpotent (all its eigenvalues are 0). Otherwise, Mat has 
eigenvalues: 

• with multiplicity 2™ - 2™, 

• 1 with multiplicity (2 n + 2(-l) n )/3, 

• j and j 2 with multiplicity (2 n — (— l) n )/3 ; 

where j = e 2l7r / 3 and n = \N/3] . 
Proof. We start from the identity 

d 

G c {t) = d + Y, Z c n {2K N)t k = tr(l - m N y l = rr;. 

k>l i=l 1 

This allows us to read off the eigenvalues directly from the generating function of the 
numbers Z^(2k, N), given by Theorem [71 

When N =3 1, G c (t) = d and so all the eigenvalues of Mat are zero. Otherwise, 

G c {t) = d + 2 n -^ + (-l) n * 



1-i 3 v ; 1-t 3 ' 
and the result follows by a partial fraction expansion. 



1 When A*' is even, the matrix Pat also describes the independent sets of the parallelogram 
V(N/2,2), as illustrated in Figure [TT] Accordingly, its powers will be used later to count inde- 
pendent sets on the parallelograms V(N/2, •). 
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In the case N =3 1, the above theorem gives an unexpected strengthening of 
Theorem This observation was communicated to us by Alan Sokal, merci a lui ! 

Corollary 10. Let N =3 1. For M > 2 1+ ^ N / 2 \ the alternating number 
Zti(M, N; C, D) of independent sets on the rectangle TZ(M, N) having prescribed con- 
ditions C and D on the two extreme diagonals of slope 1 is 0, for all configurations 
C and D. 

Proof. First assume that M = 2m + 1, so that m > 2^ N ' 2 K By the above theorem 
and the Cayley-Hamilton theorem, the mth power of M/v vanishes. Thus by j3]), the 
series Gc,D(t) is a polynomial in t of degree at most 2^ N ' 2 \ — 1. But the coefficient 
of t m in this series is precisely Zn(2m + 1, N; C, D). 

Similarly, if M = 2m + 2 with m > 2^ N / 2 \ the number Z n (M, N;C, D) is the 
entry (C,D) in the matrix M^P^, which vanishes. a 

As observed at the end of Section EJ the cylindric shape 1Z c (2k, N) can also be 
obtained by wrapping the parallelogram V(k + 1,N) on a cylinder, identifying the 
top and bottom (horizontal) layers. Consequently, the results of Theorem [7] also 
give the spectrum of another transfer matrix, denoted L^r, which describes how to 
construct the shapes V(-,N) layer by layer (Figure ITTT) . The rows and columns of 
are indexed by independent sets of the segment V(1,N). Thus the size of Ljy is 
the Fibonacci number -F/v+i, with Fo = F\ = 1 and F/v+i = Fn + Fn-i, and if C 
and D are independent sets of V(1,N), 

^N{C,D) - I q otherwise. (8) 

The generalities of Section HO imply that the spectrum (fii, . . . , fid) of Lat satisfies 

Fn+i 

F N+l +Y,Z c n {2k,N)t k = Y, T— -f 

k>l i=l ^ l 

Comparing with ([?]) shows that the spectra of L^r and Mat coincide, apart from the 
multiplicity of the null eigenvalue. 

Theorem 11. Let N > 1. The transfer matrix L^r defined by ((8|) has size Fjv+i- If 
N =3 1, then L^r is nilpotent. Otherwise, L^r has eigenvalues: 

• with multiplicity -F/v+i — 2 n , 

• 1 with multiplicity (2 n + 2(-l) n )/3 ; 

• j and j 2 with multiplicity (2 n - (— l) n )/3, 

where j = e 2j7r / 3 and n = \N/3] . 

As for the matrix Mat, the nilpotent case N =3 1 gives the following corollary, 
which has to be compared to Theorem El 

Corollary 12. Let N =3 1. Then, for K > -F/v+i, the alternating number 
Zj>(K, N;C, D) of independent sets on the parallelogram V(K,N) having prescribed 
conditions C and D on the top and bottom row is 0, for all configurations C and D. 

Proof. The number Z-p{K, N; C, D) is the entry (C, D) in the (K — l)th power of the 
transfer matrix Lat. But this power vanishes by the Cayley-Hamilton theorem. a 
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6.3. Partial results on two other spectra 

We focus in this section on two transfer matrices that generate the usual cylinder 
C(K, N) := {0, 1, ... , K-l}xZ/NZ. This cylinder can be obtained by identifying the 
diagonal borders of the parallelogram V(K, N + 1). Alternatively, it can be obtained 
by wrapping the ordinary K x (N + 1) rectangle {0, 1, . . . , K — 1} x {0, 1, . . . , iV} on 
a cylinder, identifying the vertices (i,0) and (i,N). Underlying the first construction 
are the matrices P2K defined at the beginning of Section l6~2l (see J5]) and Figure [TT1), 
Underlying the second construction is the transfer matrix Ok that describes how 
to construct the ordinary rectangles of width K. This matrix has size Fjc+i, the 
(K + l)st Fibonacci number, and its rows and columns are indexed by independent 
sets of the K point segment. If C and D are two of these independent sets, 

o^>={r Di — 01 (9) 

The similarity with the definition ([8]) of the matrix is striking, but the spectrum 
of Ok is definitely more complex than that of L^r. It is conjectured in [3] that all 
the eigenvalues of Ox are roots of unity. 

Let Zc(K,N) denote the alternating number of independent sets on the cylinder 
C(K,N). From the generalities of Section [67TT we have: 

J] Z C (K, N)t N = tr(l - tP 2K y 1 - 2 K = tr(l - tO^)" 1 - F K+1 . 

JV>1 

That is, the spectra of the matrices P2K and coincide, apart from the multiplicity 
of the null eigenvalue. 

Alas, we do not know what the numbers Zq(K,N) are. However, remember from 
Section EH] that enumerative results on configurations with open boundary conditions 
provide partial informations on the spectrum of the transfer matrix. Here, we exploit 
the results of Section [5] on parallelograms to obtain some information on the spectrum 
of P 2 K (and thus of O k ). 

For all C, D C {1, . . . , K}, let Z-p(K, N; C, D) be the alternating number of inde- 
pendent sets of the parallelogram V(K, N) having borders C and D respectively on 
the leftmost and rightmost diagonal. Then j3]) gives: 

G c , D (t) = ^-^2Kr 1 (C,D) = 6 c ,D + (-i) lCl+lDl (-l)'V 

I=(C,Ci,...,C n -i,D),n>l 

= S c ,d + H) |C|+|D| Z ^ K i N + 1;C, D)t N 

N>1 



(10) 

Note that in the first formula, the weig ht (-i)l c l+l D l(-l)l 7 l results in a weight i for 
each vertex of the extreme diagonals, as it should. Since the numbers Z-p(K,N + 
2; 0,0) coincide with the numbers Z-p(K,N) given in Theorem El this allows us to 
find some eigenvalues of P2K and Ok- We indicate in the next section how our pivot 
principle can be extended to count independent sets with prescribed borders, so as 
to determine more series Gc,d and thus more eigenvalues of the matrix 



Proposition 13. The transfer matrix Ok defined by ((H) satisfies the following prop- 
erties. 

• If K =3 1, then e i7r / 3 and e -Mr / 3 are eigenvalues o/< 
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• If K =3 2, then all the 2Kth roots of unity, except maybe — 1, are eigenvalues 
of Ok- 

• If K =3 0, then all the (2K + 2)th roots of unity, except maybe — 1 and, if K 
is odd, ±i, are eigenvalues of Ok- 

Proof. Instantiating ([TO]) to C = D = gives 

G 0)0 (t) = (1 - tP 2 ^)~ 1 (0,0) = 1 + E Zv{K,N)t N+ \ 

with Zp(K,0) = 1. Recall that the reciprocals of the poles of this series are eigen- 
values of Ok and ^2K- The numbers Zj>(K,N) are given in Theorem If K =3 1, 

G , = 1 + 1 + E(-!) n + = 1 _ t \ t 2 - (11) 

n>l 

If # =3 2, 

(2^-l)/3 2 ^ +2 

G W = 1 + E E ^^ ^^ Ma • 

Note that (1 — t 2K+2 ) /(l — t 3 ) is a polynomial, and that the only root this polynomial 
shares with 1 — t 2K is —1. 
Finally, if K =3 0, 

2K/Z-1 

g = 1 + e t 2qi - K+i)+i + e t 2qi ~ K+i)+2K+2 + E E t 2q ( K+i)+3p+3 

q>0 q>0 q>0 p=0 

1 ( ,l-t 2K 

l + t + t 3 



~ 1 _ t 2K+2 y l_ t 3 

Note that (1 + t + t 3 (1 - t 2K )/(l - t 3 )) is a polynomial. The roots it shares with 
(1 - t 2K+2 ) are -1, and, if K is odd, ±i. u 



7. Final comments and perspectives 

7.1. Other quadrangles 

A natural generalization of the rectangles and parallelograms studied in Sections [3] 
and [5] is the subgraph Q(M,N) of the square lattice induced by the points (x, y) 
satisfying 

ay < x < ay + M — 1 and — bx < y < —bx + N — 1, 

for given values of a and b. We have solved above the cases (a,b) = (1,1) and 
(a, b) = (—1,0). In particular, we have proved that in both cases, the alternating 
number of independent sets is always or ±1. What about other values of a and 
6? The case (a, b) = (0,0), which describes ordinary rectangles, shows that the 
simplicity of our results cannot be extended to all pairs (a, b). Indeed, even though 
the eigenvalues of the transfer matrix are conjectured to be roots of unity, the 
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alternating number Z{K, N) of independent sets on a K x N rectangle does not show 
any obvious pattern. For instance, for K = 4, 

1 + t 4 
(l_f2)( 1+t 3) 

1 + t 2 - t 3 + 2t 4 - t 5 + 3t 6 - 2t 7 + 3t 8 - 3t 9 + 4t 10 + O (t u ) . 

In particular, Z(4,N) ~ (— 1) N/3 as N goes to infinity. 

In contrast, recall that it is conjectured that for an ordinary rectangle with cyclic 
boundary conditions, the alternating number Zq(K, N) is 1 or —2 when N is odd [8]. 



"£z(4,N)t N = 

N>0 



(0,0) 



Figure 12. The graph £(14, 17) obtained for a = b = 2. 



Still, the simplicity of the results obtained for rectangles and parallelograms ex- 
tends to other quadrangles. For instance, if a = b = 2, a pivot rule similar to the 
one used in the proof of Theorem O (the tentative pivot lies as high as possible on 
the leftmost line of slope —2) produces the following results, where the alternating 
number of independent sets is denoted Zg(M, N). The proof is left to the reader. 

Theorem 14. Let M,N > 1. Denote m = [M/5] and n = [N/5]. 

• // M = 5 and N ^ 3, or if M = 5 1 ; or if N = 5 1,2, then E(Q(M,N)) is 
contractible and Zg(M,N) = 0. 

• Otherwise, T,(Q(M, N)) is homotopy equivalent to a sphere of dimension mn— 
1, andZ g (M,N) = {-l) mn . 

It would be worth investigating which values of a and b produce similar results. 



7.2. Transfer matrices 

We have determined in Proposition [13] some of the eigenvalues of the "hard" trans- 
fer matrix Ok- Recall that all its eigenvalues are conjectured to be roots of unity. 
We give below the value of the characteristic polynomial P{K) of the matrix O(K), 
for 1 < K < 10, and split this polynomial into the part that is explained by Propo- 
sition [13j and the (bigger and bigger) part that is left unexplained. These data have 
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been obtained with the help of Maple 
1 + t 3 



P(l 
P(2 
P(3 
P(4 
P(5 
P(6 
P(7 
P(8 
P(9 

P(10 



1 + 
1-t 4 

1 + t 
1 



1 + t 3 



1 + t 
1-t 10 



11 



1 + f. 
1-t 

1+t 3 

1 + t 
1 _+16 



1. 
1. 

(l-t 2 )(l 

(1+* 4 ), 

4\2 



(1-t 

(1 +t 4 )(l ~ t 12 ) (l - : 
1 + t 2 



1-181 



1 + t 
1 - 



(l-t 2 )(l-t 4 ) 2 (l+t 8 )(l-t 



22 i 



20 



1 + t 4 )(l - t 14 )(l + t 10 )(l - t 2U )(l - t 



,20 \ 



,26^ 



(l + t)(l + t 2 ) 
1+t 3 



(1 



1-t 2 

t 4 ) 2 (l-t 18 ) 2 (l-t 24 ) 3 (l 



t 



1 + t 



1 + t 4 



By Section 16.11 we know that each missing factor occurs in at least one of the se- 
ries Gc,c(t) counting independent sets of the parallelogram V(K,-) with prescribed 
border C on extreme diagonals. Conversely, any series Gc,z)(t) may provide some of 
these missing factors (see j3])). Hence the following question: can our pivot approach 
be recycled to compute some of these series, and do we obtain new eigenvalues in 
this way? 

For C,DC [1, Kj, denote by Z-p(K, N; C, D) the alternating number of indepen- 
dent sets of V(K,N) having border conditions C and D, respectively, on the first 
(last) diagonal. (This notation was already introduced in Section 16.31 ) Recall in 
particular the connection (flu]) between these numbers and the series Gc,D(t)- 




Figure 13. The pivot rule applied to V(4,N), with first diagonal 
{2,3} and second diagonal {4}. 

Take K = 4 and C = {2, 3}. If the configuration is C on the first diagonal, then, 
in the second diagonal, only the vertex labeled 4 may belong to an independent set. 
This gives 

Zp(4,Af;{2,3}, J D) = Z P (4, N - 1; 0, D) + Z P (A, N - 1; {4}, D) 

= Z V {4,N- 1;0,D) + Z V (4,N -4;{2,3},D). K ^> 
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The second identity is obtained by applying the pivot rule of Section to the inde- 
pendent sets counted by Zp(A, N — 1; {4}, D) (Figure [TBI. The above identity is valid 
for N > 7. We first specialize it to D = 0, and work out what happens for small 
values of N. Upon summing over N, we obtain 

= ~YZTA ( G U(t) ~t) = - {1+t){1 t _ t + t2 y 

(We have used (fTTj) for the value of G^g.) We then specialize (fT2]) to D = C = {2, 3}. 
After working out what happens for small values of N, we obtain 

r M = l-tG^dt) l + t 2 + t 3 

C ' C{ ' 1-t 4 (l + t)(l -* 4 )(1 -t + t 2 ) 1 

which now explains the missing factors (1 + i)(l — i 4 ) in P(4). 

For if = 5, we obtain similarly the missing factor 1+i 4 by considering the numbers 
Zj>(K, N; C, D), with C = {3, 4}. Indeed, after a discussion about the vertex labeled 
5 in the second diagonal and a few applications of the pivot rule, one finds 

Z V (K, N; {3, 4}, D) = Z V (K, N - 4; 0, D) - Z V (K, N - 4; {3, 4}, D), 

from which we derive 

Gc,c(t) 



(l-i 5 )(l + i 4 )' 
This gives the missing factor (1 + i 4 ). 

It would be interesting to know how far one can go with this approach. That is, can 
we determine all series Gc,c in this way? This would allow us to count independent 
sets on the ordinary cylinder, and thus to find all eigenvalues of the matrix O^. 

Note that the largest cyclotomic factor that occurs in the polynomial P{K) seems 
to be $4^-10- 

Remark added to the paper (07/03/2007): Sonja Cukic and Alexander En- 
gstrom recently pointed out to us that the following lemma [2j Lemma 2.4] can be 
used to derive some of our topological results. 

Lemma 15. Let v,w be vertices in a graph G. If N(v) C N(w) then S(G) collapses 
onto S(G — w). 

Indeed, one can match all independent sets containing w by adding or removing 
the vertex v. 

For example, for the rectangles of Theorem [6j one can eliminate the vertices in 
every third diagonal by a repeated application of Lemma [TU This leaves a graph 
formed of several paths. One can then eliminate every third vertex in each of these 
paths. The remaining graph is a disjoint union of edges, and, if M =3 1 or N =3 1, 
also isolated vertices. Thus Y,(1Z(M, N)) collapses onto an octahedral sphere (an 
mn-fold join of two points) or to the join of an octahedral sphere with a simplex, 
respectively. Theorem [6] follows. 
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